NPAC-12-01 



Complex Scalar Singlet Dark Matter: Vacuum Stability and 

Phenomenology 

Matthew Gonderinger,^'!^] Hyungjun Lim,^'[^ and Michael J. Ramsey-Musolf^' ^'[^] 

^Department of Physics, University of Wisconsin-Madison 
Madison, WI 53706, USA 

'^Kellogg Radiation Laboratory, California Institute of Technology 
Pasadena, CA, 91125, USA 

We analyze one-loop vacuum stability, perturbativity, and phenomenological constraints on a 
complex singlet extension of the Standard Model (SM) scalar sector containing a scalar dark matter 
candidate. We study vacuum stability considerations using a gauge-invariant approach and com- 
pare with the conventional gauge-dependent procedure. We show that, if new physics exists at the 
TeV scale, the vacuum stability analysis and experimental constraints from the dark matter sector, 
electroweak precision data, and LEP allow both a Higgs-like scalar in the mass range allowed by the 
latest results from CMS and ATLAS and a lighter singlet-like scalar with weak couplings to SM par- 
ticles. If instead no new physics appears until higher energy scales, there may be significant tension 
between the vacuum stability analysis and phenomenological constraints (in particular electroweak 
precision data) to the extent that the complex singlet extension with light Higgs and singlet masses 
would be ruled out. We comment on the possible implications of a scalar with ~ 125 GeV mass and 
future ATLAS invisible decay searches. 



I. INTRODUCTION 

The Standard Model (SM) of particle physics is 
known to be an incomplete theory in part because of 
its inability to explain phenomena such as dark matter 
and the baryon asymmetry of the universe. Among the 
many models vying to supplant the SM, scalar exten- 
sions of the SM are among the simplest. A gauge sin- 
glet real scalar extension has been studied as a potential 
dark matter candidate (see [rHTO]). for its impact on the 
electroweak phase transition (EWPT) [TTHH], and for 
its role in Higgs collider phenomenology [111 [15]. The 
collider phenomenology and dark matter prospects of a 
complex scalar gauge singlet [121 HZ] and a real scalar 
SU{2)l triplet have been studied as well. Exten- 

sions involving four or more new degrees of freedom, such 
as the 2 Higgs doublet model, have been widely analyzed 
over the years. Indeed, such scalar extensions of the SM 
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have resulted in a prolific field of study. 

Important assessments of the theoretical self- 
consistency of scalar extensions are the vacuum stabil- 
ity of the rcnormalization group (RG) improved one-loop 
effective potential and perturbativity of the scalar cou- 
plings. Within the SM, vacuum stability and perturba- 
tivity have resulted in theoretical bounds on the Higgs 
mass (see [52] for a comprehensive review; also, see [23]- 
[25]). In the real singlet extension of the SM, vacuum 
stability and perturbativity again place bounds on the 
Higgs mass but also constrain the singlet in dark matter 
and infiation scenarios ^29-31|. All of these results are 
dependent upon the cutoff scale of the theory, A. This 
is the scale of new physics, the scale above which new 
massive degrees of freedom can no longer be integrated 
out of the theory and become relevant for interactions 
and the effective potential. It is widely anticipated that 
new physics, in particular supersymmetry (SUSY), will 
appear at the TeV scale. The requirement of vacuum 
stability and perturbativity up to A = 1 TeV is thus a 
minimal requirement of the scalar extensions of the SM 
and result in the weakest theoretical constraints on the 
models. However, given the lack of signatures of SUSY 
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or other new physics in early LHC data (see for exam- 
ple [32H35] ). it is possible the scale at which new physics 
and new massive degrees of freedom become relevant lies 
beyond the TeV scale. As the cutoff scale increases, vac- 
uum stability and pcrturbativity can impose increasingly 
significant constraints on the scalar extensions. 

In this work, we study this issue of vacuum stability 
and perturbativity — particularly for higher cutoff scales 
— for the complex scalar singlet extension of the SM, re- 
ferred to as "the CccSM". With an appropriate set of 
symmetries, this model yields both a viable dark matter 
candidate {A) as well as two real neutral scalars h' and 
S' that are mixtures of the SM Higgs boson and the real 
part of the complex singlet. We discuss in detail the re- 
quirement of vacuum stability of the effective potential, 
i.e., that the clectroweak minimum of VcS be deeper than 
any other minimum. Generally, however, there exists at 
least one deeper minimum at large values of the scalar 
field (p due to top quark loop contributions. The con- 
ventional stability requirement, then, is to restrict the 
effective theory to energy scales below the value of jp 
for which Vcs{ip) falls below the electroweak minimunrj 
One then identifies the maximum stability scale A with 
this maximum value of ip. This criterion, however, is 
gauge-dependent since only the value of the potential at 
its extrema is gauge-invariant . whereas the field itself 
remains gauge-dependent. Consequently, identifying the 
cutoff A with a value of ip is not physically meaningful. 
As an alternative, we will use an analysis of the RG evo- 
lution of the quartic couplings of the CxSM since the 
effective potential is dominated by terms quartic in the 
fields. By restricting the running of these parameters to 
energy scales below an appropriately chosen value of A, 
we guarantee in a gauge-invariant way that the effective 
potential is bounded below and that the EW minimum 
is stable. 

In addition to stability and perturbativity considera- 
tions, we will also apply various phenomenological con- 
straints in our analysis of the CxSM: results for elec- 
troweak precision observables (EWPO), dark matter relic 
density and direct detection measurements, and limits 
from LEP. We also study scenarios that may lead to a 
strong, first order electroweak phase transition (EWPT) 
as is needed for electroweak baryogenesis that may lead 
to relic gravity waves. We find that, should new physics 
exist at the TeV scale, the CxSM has regions of param- 
eter space which satisfy all constraints and favor a rela- 
tively light and weakly coupled singlet-like scalar in ad- 
dition to a Higgs-like scalar in the current mass range 
allowed by searches at ATLAS [S^ and CMS [38]. Con- 
versely, if new physics does not appear until higher en- 
ergy scales well above a TeV, the vacuum stability con- 



Alternately, if the electroweak minimum is metastable with a 
lifetime longer than the age of the Universe, then the stability 
radius in field space can be increased. 



siderations are in significant tension with experimental 
constraints, particularly EWPO data. Rather generally, 
the CxSM can be ruled out should new physics fail to 
appear just below the grand unification (GUT) scale, 
Mqvt — 10^^ GeV. These conclusions hold for both a 
relatively light dark matter mass (as indicated by the Co- 
GeNT [SHU], DAMA/LIBRA gHH], and CRESST-II 
[15] collaborations) and a heavier dark matter mass. Fur- 
thermore, it is possible for the scalars h' and S' to decay 
to dark matter; in fact, in scenarios where both the dark 
matter and one of the scalar eigenstates are light and the 
scale of new physics is roughly a TeV, the branching frac- 
tion to dark matter is sufficiently large that the ATLAS 
detector could be sensitive to these invisible decays. 

Our discussion of these issues is organized as follows. 
We begin with an introduction to the CxSM in section[lT| 
In section |III[ we describe the requirements of vacuum 
stability and perturbativity in detail, and discuss in de- 
tail the impact of gauge dependence on the traditional 
vacuum stability analysis. We then present our analysis 
of the RG evolution as a gauge-independent substitute. 
Section |IV| introduces phenomenological constraints on 
the CxSM from the EWPT, EWPO, dark matter rehc 
density and direct detection measurements, and collider 
physics at LEP and the LHC. We present our results in 
section [V] Section IVll contains our conclusions. 



II. COMPLEX SINGLET MODEL 

II. 1. Tree Level Potential 

In the CxSM, the SM is supplemented by the addi- 
tion of a single complex scalar degree of freedom that 
transforms trivially under the SM gauge groups. Thus, 
the only renormalizable tree-level interactions between 
the complex singlet, §, and the SM occur in the scalar 
potential of eqn. ([T]) — the singlet couples to the SM 
fermions and gauge bosons only through the Higgs, H 
(sometimes referred to as the "Higgs portal" |;46]). 

V{H,S) = ^m^H^H+^ {H^nf 

-i-|ijti/|§|V||§r + ||§r 
+ Q e'^"S2 + |ai| e^^-^S + c.c.^ (I) 

In the absence of the bi and ai terms, V {H, §) obeys 
a global C/(l) symmetry: S e'"S. By breaking this 
symmetry both spontaneously and softly (through the 
last two terms), we obtain a cold dark matter can- 
didate. When the singlet gets a vacuum expectation 
value (vev), (§) = x/\/2 (the Higgs has its usual vev, 

(H) = {0,v/V2)^ where v = 246 GeV), the global ;7(I) 
symmetry is spontaneously broken, the real part of the 
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singlet mixes with the SM Higgs, and the imaginary part 
of the singlet becomes a massless Goldstone boson. To 
give mass to the imaginary part of the singlet so that it 
can potentially fill the role of a stable cold dark matter 
candidate, we include the explicit J7(l)-breaking terms 
proportional to bi and oi. 

Note that for ai = the potential retains a Z2 sym- 
metry associated with the components of S. Since spon- 
taneously broken discrete symmetries create issues with 
cosmological domain walls [TTHIS] we also introduce an 
explicit Z2-breaking term proportional to oi. These ad- 
ditional terms are chosen so that the potential retains 
a Z2 symmetry for Im(§), thereby ensuring stability of 
the dark matter particle. Moreover, these operators close 
under renormalization. The phase (f>ai can be absorbed 
in a redefinition of S and 4>bi, and we choose (f>bi — ti" 
to avoid mixing between the real and complex compo- 
nents of S [I7|. Then, expanding S = (S" + iA) /V2 antQ 
H = h/\/2 gives the tree level potential 



Vo {h, S,A)^'^h^+ + |/^2 (^2 ^ ^2) 

-t- J (62 - bi) 52 + i (62 + bi)A^ - V2aiS 



'S'A' + ^{S^ + A^) . (2) 



Requiring that the potential in eqn. ^ have a mini- 
mum at (H) = h/V2 = v/V2a.nd (§) = S + iA = x + i-0 
gives the following set of minimization conditions: 



dVo 
dh 



0, 



dS 



dA 



(3) 



where all derivatives are evaluated at (/i, 5, A) = (v, x, 0). 
(Note that other solutions to the minimization equations 
may exist; however, our vacuum stability analysis de- 
scribed in section ITTll verifies none of these other critical 
points is a global minimum given values for all of the 
parameters.) These minimization conditions allow the 
Higgs vev V and the singlet vev x to replace and 62 
as parameters in the CxSM according to eqn. Q. 



2 _ 1 \ 2 

m = Av ' 

2 

62 = fei + 2V2 



fli 



d2x'^ 



(4) 



At the minimum, the mass (second derivative) matrix is 



then given by eqn. ([s]). 
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^Is 


"^Ia 


^Is 


m| 


rn^A 






m\ 



\52XV 

\52xv ^d2X^ + \/2ai/x 




bi + ^/2ai/x 



(5) 



We choose the ^7(1) and Z2 symmetry breaking parame- 
ter ai such that ai <^ x — i.e., we take ai = 10^^ GeV^ 
and a; > 10 GeV. This choice serves two purposes: first, 
it simplifies the model by reducing by one the number 
of unknown parameters that must be varied; second, it 
ensures that the minimum at {h,S,A) = {v,x,0) is the 
global minimum of the potential, as we will discuss in sec- 
tion |TTT] With this choice for oi, the dark matter mass 
is given by ~ y/bi- Meanwhile, the non-zero entry 
for mlg induces mixing between the SM Higgs and the 
real component of the singlet. The resulting mass eigen- 
states, which we denote h' and S", have masses given by 
the eigenvalues of M, the 2x2 upper left quadrant of 
eqn. These eigenvalues are 



Tr (M) ± J (Tr (M))^ - 4Det (M) 



(6) 



with m-|_ > m_. In order for these masses to be positive 
real numbers, the condition Det (M) > must hold. In 
the limit of small ai , this condition simplifies to (5| < \d2 . 
The eigenstates h',S' are written in terms of h and S 
according to eqn. Q. 



(7) 



The eigenstates h' , S' couple to the fermions and gauge 
bosons via SM Higgs couplings reduced by a factor of 
cos 4>, — sin (f), respectively. The mixing angle (p is given 
at tree level by 



h' 




COS (p 


sin (j) 




h 


S' 




— sin 


COS0 




S 



tan 2(/) = 



(8) 



We take the mixing angle to be —tt/4 < (j> < tt/A so that 
h' is always the "Higgs-like" eigenstate and S' is always 
"singlet-like'l^ Which eigenstate is heavier will depend 
on our choice of parameters. When choosing 62 < 



^ We always choose the minimum of the potential so that the neu- 
tral real component of the Higgs doublet has a non-zero vev and 
the other components, the would-be Goldstone bosons, have zero 
vev. 



In the literature, the mass eigenstates are often denoted as hi 
and h2. We use a different notation to emphasize that one state 
is always "Higgs-like" and the other is always "singlet-like" . 
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(which allows for a strongly first order EWPT while be- 
ing consistent with LEP bounds, as discussed later in sec- 
tion IV.l ), eqns. ([5])-([8]) imply that h' will be the heavier 
eigenstate with < for relatively large A and relatively 
small c?2, X whereas S' will be the heavier eigenstate with 
> for relatively small A and relatively large ^2,2;. 



II. 2. One-Loop Potential 

For our vacuum stability analysis, we use the full 
Coleman- Weinberg one-loop effective potential at zero 
temperature with one-loop renormalization group (RG) 
running parameters. 



Kff {h, S, A) = Vo (/i, S, A) + Vi {h, S, A) 



(9) 



Vo {h, S, A) is given in eqn. ([2]), where all fields, couplings, 
and masses are replaced by their RG running counter- 
parts. The one-loop contribution, calculated in the Lan- 
dau gauge and renormalizcd in the MS scheme, is given 
by 



Viih,S,A) = 
1 



647r2 



n,Tr log 



7^ 



(10) 



The sum i runs over scalars, fermions, and gauge bosons. 
The field-dependent mass matrices Mf , the number of 
associated degrees of freedom n^, and the numerical con- 
stants Ci are given in appendix|X] /i is the 't Hooft renor- 
malization scale. As discussed in the next paragraph, 
the effective potential is renormalization scale indepen- 
dent to one-loop order; any residual scale dependence 
is higher order. To remove this residual scale depen- 
dence, we would like to choose ^ to minimize the logs in 
eqn. ( 10 ). However, no single choice for /i will simultane- 



ously minimize all of the logs, and so we make the simple 
choice = + S'^ + A^ . 

The RG equations for the fields, couplings, and 
masses are determined by requiring scale invariance of 
the effective potential to one-loop order: the scale de- 
pendence implicit in the parameters of Vo cancels the 
explicit scale dependence in Vi, i.e.. 



a/i dfi Ofi 



(11) 



Applying this condition to the CxSM effective potential 
gives a series of equations to be solved for the (5 and 
7 (anomalous dimension) functions that determine the 
running of the fields, couplings, and masses. The j3 and 
7 functions are given in appendix |Xj For convenience, 
we take /i = Mz as the input scale for all our running 
parameters. 

In analogy with the tree level potential, we apply the 
minimization conditions to the one-loop effective poten- 



tial: requiring that the minimum of the effective po- 
tential occur at (/i, S, A) = {v, x, 0) fixes the bound- 
ary conditions for the running mass parameters (Mz) 
and &2 (Mz). Furthermore, we obtain the masses of the 
scalars by diagonalizing the matrix of second derivatives 
of the effective potential evaluated at the minimum. The 
dark matter A is protected by a Z2 symmetry so that it is 
stable and does not mix with h and S at the minimum of 
the potential even upon inclusion of the one-loop correc- 
tions. (The necessity of the Z2 symmetry to ensure the 
dark matter cannot decay is a generic feature of these 
simple scalar extensions of the SM. The real scalar sin- 
glet dark matter extension of the SM is referred to as 
the Z2a::SM for this reason.) The mass eigenstates h' , S' 
are defined in terms of h and S as in eqn. ([7| using the 
one- loop value of the mixing angle 0. 



III. VACUUM STABILITY ANALYSIS 



III.l. The Vacuum Stability Analysis 



The requirement of absolute vacuum stability is 
equivalent to requiring that the electroweak (EW) zero- 
temperature minimum of the effective potential be a 
global minimum over the energy range for which the SM 
is valid. The common practice for vacuum stability anal- 
yses in the literature begins, as described above, with 
the RG improved effective potential (generically VcS {'Pi)) 
and choice of the renormalization scale fx^ ^ 0^ = fifi 
to minimize logarithms in the one-loop potential. Then, 
the maximum radius in field space, ^max, is found ac- 
cording to the requirement of absolute vacuum stability: 



EW) 



(12) 



where 0ew gives the values of the fields at the elec- 
troweak minimum — (/i, S', A) = (w, x, 0) in the CxSM. 
This maximum radius in field space is identified with the 
cutoff scale of the effective theory, A. It's presumed that 
A = 'Pmax is the scale at which new physics is required to 
alter the shape of the potential so the electroweak mini- 
mum remains a global minimum. 

The requirement of absolute vacuum stability can be 
relaxed to the case of metastability for which the EW 
minimum may not be a global minimum, but the tunnel- 
ing probability from the EW minimum to the true global 
minimum is sufficiently small (the lifetime of the elec- 
troweak vacuum is greater than the present age of the 
universe). In the real scalar singlet extension of the SM, 
the authors of [50] showed that the vacuum metastability 
requirement is indeed less restrictive of the model param- 
eter space than the absolute vacuum stability analysis in 
|29) . To obtain more conservative bounds, and for sim- 
plicity, we focus on the absolute stability scenario. 
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III. 2. Vacuum Stability for Scalar Extensions of 
the SM 



There are two primary considerations for vacuum sta- 
bility in the CxSM. The first is the possibiUty of a Z2 
symmetry breaking minimum at tree level; the second is 
the set of constraints on the quartic couplings and their 
RG evolution. 



ni.2.1. Z2 Symmetry Breaking Minimum 

In the Z2XSM, the singlet mass depends on its 
quadratic mass parameter and its coupling to the Higgs. 
In [W, it was shown that for dark matter masses in the 
range 10-100 GeV there is a tension between having a suf- 
ficiently large (positive) Higgs-singlet coupling to avoid 
oversaturating the dark matter relic density and main- 
taining a stable EW minimum of the potential. Obtain- 
ing light scalar singlet dark matter that saturates the 
relic density can require a negative mass-squared param- 
eter in the potential, leading to a minimum along the sin- 
glet axis of the potential for which {S) ^ and {H) — 0; 
thus the Z2 symmetry is broken and the dark matter 
is not a stable particle. In the present analysis of the 
Ca;SM, the dark matter mass depends on the linear pa- 
rameter ai and the quadratic parameter bi (plus small 
loop corrections). We have chosen ai to be small, and so 
to obtain a positive dark matter mass we unambiguously 
choose 61 to be positive. Thus there is no dangerous 
symmetry-breaking minimum along the Im [S]-axis and 
we do not have a tension between the dark matter mass 
and vacuum stability as in the real scalar singlet model. 



III. 2. 2. Limits on the Quarttc Couplings 



In the CxSM, the stability of the tree level potential 
minimum is guaranteed simply by requiring thai|^ 



S2 < Ad2 , 
A > , 
d2 > . 



(13) 



The first condition is necessary for obtaining positive 
mass-squared eigenvalues for the mixing between the 
Higgs and real component of the singlet at the minimum 
of the potential, {h,S,A) = {v,x,0). This is of course 
equivalent to the second derivative test to ensure that 
the critical point {h,S,A) — (w,x, 0) is actually a min- 
imum. The second two conditions are required for the 
potential to be bounded below in all scenarios; the first 



condition is also required for the potential to be bounded 
below when ^2 < 0. 

Going beyond tree level with the one-loop potential 
and the one-loop RGEs affects the stability of the poten- 
tial in two ways. 
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As in the SM and Z2XSM, a global minimum for 
h ^ V can arise due to the running of the Higgs 
quartic coupling A. The large top Yukawa coupling, 
Ut, causes A to evolve to negative values for large 11 
when A (Mz) is sufficiently small, as seen from the 
/3-function in eqn. (14). 



I3\ = M 



(14) 



Along the ft,-axis of the potential, the Higgs self- 
coupling dominates, viz.: 



ft. > 



(15) 



Setting fj, — h to minimize large logs in the effec- 
tive potential thus combines these two effects: the 
potential can develop a very deep global minimum 
(/i > u, = 0, A = 0), high above the EW scaleQ 

In the Z2XSM, the Higgs-singlet coupling 62 has a 
positive contribution to the /3-function of A irre- 
spective of the sign of S2 ■ As was shown in [29] , the 
contribution of 82 to the running of A decreases the 
theoretical lower bounds on the Higgs mass from 
vacuum stability. In the CxSM, the running of A is 
again tempered by the Higgs-singlet coupling ^2, so 
a larger value of S2 may push this deep minimum 
above the cutoff scale A. 

The second effect is one specific to the choice 62 < 
when one-loop corrections are included. At tree 
level, the requirement 5| < Ac?2 is sufficient to pre- 
vent a runaway direction in the potential between 
the h-, S-, and A-a.xes when 62 < 0. However, as 



the /3-function for 62 shows (eqn. (v416)), a negative 
82 {Mz) will decrease as the scale fi increases. This 
could in principle lead to a runaway direction in 
the potential for some region of field space between 
the axes. 



III. 3. Gauge Dependence 

It has been pointed out that the SM one-loop effective 
potential depends on the choice of gauge- fixing condition; 
equivalently, in the i?^ gauges, the potential depends 



* See | 17| for further discussion. 



^ This is a minimum and not simply an unbounded direction be- 
cause A does become positive again at higher scales. 
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Ma=10 GeV, x=100 GeV, ^2=0.2, A=10^3 GeV 
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FIG. 1. A plot of (5i (Mz) vs. A (Mz) d2 (Mz). For all points, 
Ma = 10 GeV, x = 100 GeV, and d2 (Mz) = 0.2. The tree 
level vacuum stability requirement, S2 < Ad2, is indicated 
with the solid line. All points satisfy the effective potential 
vacuum stability requirement of eqn. ( 12 1 with A = 1 TeV 
for some choice of gauge parameter ^. For gray points, ^ = 
(Landau gauge); for blue points, ^ = 1, and for red points 
^ — 50. (In this and subsequent figures, the point (0,0) is 
included for reference only.) 



on the gauge parameter ^ (see [STJ [52] and references 
therein). Hence, the field expectation values at a mini- 
mum, ifmin, do not correspond to a physical observable. 
It was shown in [51] that the value of the effective poten- 
tial at its extrema can be calculated in a gauge-invariant 
way through a consistent expansion in h, provided that 
the extrema have classical (tree level) analogs. The va- 
lidity of this procedure is a consequence of the Nielsen 
identity [5S] . 

The gauge dependence of the effective potential 
presents complications for the vacuum stability analysis 
despite the existence of the Nielsen identity. As described 
above, the vacuum stability analysis generally performed 
in the literature is interested in a particular radius in 
field space, ipmax, obtained from eqn. ([T2|). Because the 



potential is gauge dependent, it is possible that for one 
choice of gauge the potential may satisfy the stability re- 
quirement below ipmax, but for another choice of gauge 
the potential may become unstable: 



2 

max ' 



Vcs (</5i;Ci) > Vcs {0Ew;^i) V < , ^^^^ 



This ambiguity is dramatically demonstrated in fig. [T] for 
one particular choice of the Ca;SM parameters. The plot- 
ted points are those allowed by the vacuum stability re- 
quirement of eqn. ( 12 ) (with a 1 TeV cutoff) for three 
different gauge parameters: C = (gray), C = 1 (blue), 
and ^ = 50 (red). Thus, identifying the cutoff scale of the 
effective theory — a physical, gauge-independent number 



— with <^rnax IS problematic]^ 

There exist in the literature two methods for perform- 
ing a gauge-independent analysis of the vacuum stabil- 
ity and corresponding Higgs mass bounds: the "physical 
effective potential" in [53], and the Vilkovisky-DeWitt 
formalism in [M]. These methods have been applied to 
toy models and have derived gauge-independent results 
that reproduce to within a few percent the results of a 
traditional vacuum stability analysis done in the Landau 
gauge [Sni[5S]- To our knowledge, however, no gauge- 
independent method for analyzing the vacuum stability 
of the effective potential in the full SM, much less the 
CxSM, has been presented. 

Rather than generalizing either of the above meth- 
ods, in the current analysis we choose to make vacuum 
stability arguments based on the running of the quartic 
couplings that dominate the potential, i.e., 

+ S2 (m) {S^ + A^) + d2 ip) {S^ + A'') (17) 

for /i^ > w2[7] ti^g CxSM, vacuum stability requires 
that the tree level couplings obey eqn. (13 1. We ex- 



tend these requirements to the one-loop RG running cou- 
plings, as in eqn. ([l8|. 



(/i) < A ifi) d2 (m) 

x{^i)>o 

d2 (Ai) > 



> V /i < A 



(18) 



After evolving all the parameters of the theory in the 
effective potential up to the cutoff scale of the theory, if 
any of these conditions is violated then the potential may 
become unstable for larger scalet]^ 

• If A (/i) there will be a deep second minimum 
along the /i-axis of the potential. 

• If d2 in) 'i' there will be a "runaway direction" 
of the potential along the S- and ^-axis, i.e., the 
potential is unbounded from below. 



^ Furthermore, in the SM and its extensions such as the CxSM, 
the appearance of a global minimum for h ^ v occurs due to 
the RG running of the Higgs self-coupling A at one-loop. There 
is no classical minimum corresponding to this new global mini- 
mum appearing at one-loop, and so the perturbative h expansion 
described in 1511 yields trivial equations when evaluating the po- 
tential at this minimum in a gauge-independent way. 

^ A full gauge-invariant vacuum stability analysis of the effective 
potential in the SM and its scalar extensions is relegated to future 
work. 

* Minimizing logarithms in the one-loop potential requires the 
choice fi^ = when evaluating the potential in a gauge- 

independent fashion at the minimum 0rn.im as discussed in 
|III.2.2[ 
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If 62 (m) it A {^) d2 (/i) and ^2 < there will be 
a runaway direction somewhere between the field 
axes. 



• If (52 > 0, the requirement (5| (/i) < A {^i) 6.2 (/i) may 
be overly restrictive. Though (5| < Ad2 is necessary 
to ensure that (h,S,A) = {v,x,Q) is a minimum, 
the running of 62 at large scales will not affect the 
shape of the potential at the electroweak minimum. 
Thus the EW minimum will remain the global min- 
imum and the potential will be bounded below for 
large values of the field. 

The values of the fields where these instabilities occur is 
immaterial to our analysis; the mere fact that they occur 
because the conditions of eqn. ( 18 ) are violated implies 



that the vacuum stability requirement is not satisfied. 
Since the RG evolution of all the mass and coupling pa- 
rameters in the CxSM is gauge-independent, the scale at 
which any one of the requirements of eqn. ( 18 ) is violated 



— which we identify with the cutoff scale of the effective 
theory, A — is also gauge- independent. We emphasize 
that the constraints placed on the couplings (and hence 

are moti- 
Since we 



18f 
ity. 



the masses of the scalar fields) from eqn. ( 
vated by the requirement of vacuum stabi 
have not calculated gauge-independent tunneling proba- 
bilities for transitions to a non-EW global minimum of 
the potential, our analysis may give more conservative 
bounds than those determined by allowing the EW min- 
imum to be metastable. 



III. 4. Perturbativity 



III. 5. Analysis Procedure 



In practice, we take as inputs the boundary condi- 
tions for the running Lagrangian parameters (with the 
boundary conditions for m? and 62 fixed by the other in- 
puts, the scalar vevs v = 246 GeV and x). We then solve 
the RGEs up to the Planck scale {O (10^^) GeV) and 
determine the scalar masses and mixing angle by diag- 
onalizing the matrix of second derivatives (all evaluated 
at {h,S,A) = (w,x,0)): 



dlVcs dhdsVcs dhdAVcs 



P-Diag(A/2,,M2,,Afl) -P-i (20) 

where the matrix P is the orthogonal matrix containing 
the mixing parameterized by the angle between the SM 
Higgs and the real part of the singlet. We take the cutoff 
scale A as an output, obtained by finding the minimum 
value of /i for which either eqn. ( 18 ) is violated or one 



or more of the couplings becomes non-perturbative ac- 
cording to eqn. (19). We also scan over field space for 



the minimum radius (if such a point exists) at which the 
value of the potential is equal to its value at the elec- 
troweak minimum — V^s (h, S, A) = Ves{v,x,0) — to 
compare the gauge-independent results from the running 
couplings with the gauge-dependent results of the effec- 
tive potential in the Landau gauge. The cutoff scale is 
required to be at least one TeV. 



We also require that the couplings in the scalar poten- 
tial remain perturbative for all values of the scale fj,. The 
definition of "perturbative" is somewhat subjective. At 
one-loop order in perturbation theory, the quartic scalar 
couplings all have Landau poles as fi approaches A^p; 
minimally, the location of the Landau pole could be taken 
as the cutoff scale of the theory, A — A/,p. However, the 
couplings reach unreasonably large values well before the 
Landau pole. Two-loop analysis of the SM RGEs shows 
that Higgs quartic self-coupling A has a fixed point at 
large scales where /3a — >■ and A (^) — Xpp Fur- 
thermore, it has been shown in [57j that the SM remains 
perturbative for values of the Higgs quartic self-coupling 
A (A) in the range App/4 to App/2. A full two-loop anal- 
ysis of the Ca;SM is beyond the scope of our current 
work, so we impose the an approximate perturbativity 
constraint on the couplings in eqn. (19 1. 



S2 (m) < App/3 
A M < Xfp/S 
d2 {y) < App/3 



V Af^ < M < A 



(19) 



IV. IMPLEMENTATION OF OTHER 
CONSTRAINTS 



IV. 1. Electroweak Phase Transition 



The electroweak phase transition in a real scalar sin- 
glet extension of the SM has been studied extensively 
[TTHI3]. If the potential possesses a Z2 symmetry, as 
in the Z2a;SM, which is spontaneously broken by a non- 
zero singlet vev, the Higgs and singlet mix and the sin- 
glet cannot act as a stable dark matter candidate. Under 
these circumstances, a negative value of the Higgs-singlet 
coupling parameter, 82, is most suitable for obtaining a 
strong first-order phase transition and satisfying the LEP 
constraints on the Higgs mass and mixing angles. We 
therefore consider 82 < Q for our analysis. As mentioned 
in section [hi] this choice has the most interesting impli- 
cations for vacuum stability of the full one-loop potential 
with nmning parameters [58| . 
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IV. 2. Electroweak Precision Observables 



In the SM, measurements of electroweak precision ob- 
servables (EWPO), such as Z° pole measurements, pro- 
vide sensitivity to the Higgs mass via loop-level effects. 
In the CxSM, mixing between the SM Higgs and the real 
component of the complex singlet alters these loop-level 
effects. To determine which values of the CxSM param- 
eters best match EWPO data, we follow the procedure 
described in |12j . which we summarize here. 

The EWPO data are parameterized in terms of the 
oblique parameters S, T, and U. Experimental values of 
the oblique parameters are determined by performing a 
best fit analysis using all electroweak precision data, as 
in |12j . Alternatively, the oblique parameters can be cal- 
culated analytically as they are defined in terms of the 
self-energy corrections to the gauge boson propagators: 
n^z (p^) , ^ww {P^) , (P^) , {p^) ■ Given the di- 
rect search limit from LEP on the Higgs mass, a SM ref- 
erence value, 0° = O {Mf^^ = 114.4 GeV), can be com- 
puted for each of these oblique parameters (O = S, T, U) . 
The best-fit value determined from electroweak precision 
data for the difference between O and the SM reference 
value is defined as 



AO" = O - 



(21) 



Since the real component of the scalar singlet § mixes 
with the Higgs, the propagator corrections Hvi/vf and 
Tlzz in the CxSM, and hence the oblique parameters, 
are different from the SM results (however, H^^ and H^^ 
are unchanged because the scalars are neutral). The dif- 
ference can be written as 



AO 



fSM 



• O M^" Ms, 



O" 



(22) 



Thus, given values for the masses of the scalar eigenstates 
h' and 5' and the mixing angle (j), all extracted from the 
effective potential, the quantity AO can be computed. 
The masses and mixing angle of the CxSM are consistent 
with EWPO data if the oblique parameter differences AO 
fall within the 95% C.L. region of the experimental values 
AO*'. This is equivalent to Ax^ < 7.815, where Ax^ is 



defined in eqn. (23 1 using the correlation matrix p and 



errors a from [12j . 
Ax' ^ J2 (^O. - AO°) {apa)-/ (AO, - AO°) (23) 

id 

The analytic forms of the oblique parameters are given 
in appendix [B] 




, h', S' A ^ h'. S' 



^ h', S' A 



. h'. S' 



FIG. 2. Feynman diagrams showing processes contributing to 
the annihilation cross section of the dark matter particles, A. 



complex singlet is a stable dark matter candidate. The 
thermal relic density of the scalar A, ^Ah'^, is controlled 
in part by the annihilation cross section of the dark 
matter particles, aann (AA XX): roughly, ^ 
^/{o'annVrei) where {(TannVrei) IS the thermal average of 
the annihilation cross section times the relative velocity 
of the dark matter particles in the center-of-mass frame. 
Processes that contribute to <Jann are shown in fig. |IV.3| 
The kinematical/mass-dependent factors associated with 
the cross sections for these diagrams can be of particu- 
lar importance in determining the relic density. In the 
limit of non-relativistic dark matter where i/s ~ 2mA 
and Vrei ~2\pA\/mA, 



d 

dn 



(Tann {AA -> XX) Vrel OC 



(24) 



In the case of the four-point vertex in fig. IV. 3[ the am- 
plitude is independent of masses and momenta, so 
the annihilation channels AA — > h'h' or AA — > S'S', 
when kinematically allowed, are largest for Mh' or Mg/ <C 
Ma- For the s-channcl resonances AA — )■ /i' — >■ XX 



and AA 
(4Ml 



S' 



XX, \M\' 



{AMI 



Mi: 



or 



Mg,) ; thus the annihilation cross section of 
course becomes large for Mh' or Ms' — 2Ma- When 
the annihilation cross section becomes large, the complex 
singlet undersaturates the total dark matter relic density, 
^DAih^i for which we use the WMAP la measurement 
^DAih'^ — 0.92 — 0.118 [Snj. We use the computational 
tool micrOMEGAs [BD] to numerically calculate the relic 
density. Though we allow the CxSM to undersaturate 
the relic density, oversaturation is forbidden. 



IV.3. Dark Matter Relic Density 



IV. 4. Dark Matter Direct Detection 



As described in section [Hj we choose parameters for 
the CxSM such that the imaginary component of the 



A number of experiments have performed searches for 
dark matter scattering off atomic nuclei and have pub- 
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lished limits on the spin-independent scattering cross sec- 
tion as a function of the dark matter mass. The most 
restrictive hmits at present come from the XENONIOO 
[5T| i2] and CDMS [Ml El] experiments. In apparent 
conflict with these limits are results from the CoGeNT 
[SSSH, DAMA/LIBRA [42H44|, and CRESST-II \45\ ex- 
periments which have observed signal events correspond- 
ing to dark matter particles with Ma ~ 10 GeV|5S]. For 
our analysis, we remain impartial in the debate over these 
experimental results. In the CxSM, the scattering cross 
section of the dark matter candidate with a proton is 
calculated (see [TSIISS]) according to eqn. (25). 



CTdd 



27ri;2 (mp + Ma) 

f QAAh' COS (f) gAAS' sin ( 



Ml 



Ml 



fpu + fpd + fj 



ips 



27 



(3/g) (25) 



where 



9AAh' = {S2VCOS(j) + d2xsiii4}) /2 (26) 
gAAS' ^ id2X cos (j) - d2vsiii(j)) /2 (27) 



The proton matrix elements /, 



™p/t? = {plfnqqqlp) 



f(p) 
Jtg 



E 



f^^l , (28) 



are calculated in [57]; we take the central values 

f^^ = 0.020 f^J = 0.026 f^] = 0.118 (29) 

We consider masses and cross sections that satisfy ex- 
actly one of the direct detection experiments: cither the 
upper bound from XENONIOO or the signal regions from 
CoGeNT, or DAMA/LIBRA, or CRESST-II. We utilize 
micrOMEGAs [5H] for numerical calculation of the di- 
rect detection cross section. In comparing this calcu- 
lated scattering cross section to the limits from the cited 
experiments, in eqn. (30) we scale the cross section by 



the fraction of the total relic density constituted by the 
CxSM dark matter candidate to account for the reduced 
flux of dark matter particles in the detectors when the 
relic density is undersaturated. 



aled ^dd 



(30) 



IV. 5. LEP Mixing Angle Constraints 

Application of the LEP limits (this section) and the 
ATLAS invisibly decaying Higgs search conditions (sec- 



tion IV. 6) requires calculation of the scalar mass eigen- 



states' widths. These are given by eqn. (31) 

„2 



Ttot {h') = cos' (p ■ FsM [H* 
[+ F [h' S'S')] 

Ttot{S')^sm^<f>-TsM (H*) 
[+ F {S' ^ h'h')] 



r{h' 
r{h' 

T{S' 

r{s' - 



AA)] 
AAAA)] 
AA)] 
AAAA)] 



(31) 



In eqn. (31), Fsm (H*) is the rate of decays of the SM 



Higgs to SM final states where the Higgs is assumed to 
have a mass equivalent to that of the h' or S' eigenstate 
appropriately. We calculate the SM Higgs width as a 
function of the Higgs mass using the program HDECAY 
[69 . The decay rates in square brackets in eqn. (31 1 are 



only included when the indicated decay is kinematically 
allowed. The decay rate of the h' , S' eigenstates to pairs 
of dark matter particles is given in eqn. (32). 



r{h',S' AA) = 



\9AAh',AAS' r 

'i2iimh'S' 



'1-4- 



''DM 



(32) 



The parameters gAAh' and gAAS' a-re defined in sec- 
tion |IV.4| Decays to four dark matter particles have 
intermediate states of two (possibly off-shell) scalars. 

The LEP Working Group for Higgs Boson Searches 
has made use of the combined data from the four LEP 
experiments to constrain the mass and ZZH coupling of 
BSM Higgs- like scalars f70|. An upper bound is set on 
the quantity 



f:'2 — ( 9Wzh \ 
^ - \9Tzh) 



X Br{H ^ SM) 



(33) 



as a function of the Higgs mass. If the scalar particle 
H has only SM decays, then Br {H -> SM) = 1 and the 
limits are on the BSM-to-SM ratio of the Higgs- Z-Z cou- 
pling. In the CxSM, the ratio of the couplings is cos^ (p 
for h' and sin^ for 5". If additional scalar decays are 



kinematically allowed, then the widths in eqns. (31 )-(32) 
are used to calculate Br {H — > SM). We apply the LEP 
limits to both scalar mass eigenstates h' and S' . 



IV. 6. ATLAS Sensitivity to an Invisibly Decaying 
Higgs 



The mixing between the SM Higgs and the real com- 
ponent of the complex singlet, as well as the potential 
for one or both eigenstates to decay to an even num- 
ber of dark matter particles when kinematically allowed, 
also has implications for Higgs searches at the LHC. The 
study in [71] found that the ATLAS experiment at the 
LHC would be sensitive, via the vector boson fusion chan- 
nel, to a Higgs with a mass between 114-200 GeV and an 
invisible decay mode so long as the condition 

^2 = Br{H^ mvis) x > 60% (34) 

O'SM 
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is satisfied (for masses greater than 200 GeV, the re- 
quirement increases to > 70%). In eqn. (34), ctbsm 



and CTsM are the Beyond-the-Standard-Model and SM 
production cross sections, respectively. In the CxSM, 
cbsm CTSM where the proportionality factor is either 
cos^ (j) for the "Higgs-hke" eigenstate or sin^ (j) for the 
"singlet-like" eigenstate. The invisible decay branching 
fraction, Br (H — > invis), includes the kinematically al- 
lowed decays of the h' or S' to two or four dark matter 
particles: 





Ma = 10 GeV 


Ma = 100 GeV 


x = 10 GeV 


d2 = 0.2,0.5,0.9 


d2 = 0.2,0.5,0.9 


X = 100 GeV 


da = 0.2,0.5,0.9 


d2 = 0.2,0.5,0.9 


X = 1000 GeV 


d2 = 0.2,0.5,0.9 





TABLE I. A list of parameter scans performed, ai is chosen 
to be 10"'^ GeV^ so 6i ~ Ma- The only free parameters are 
A and 62- A number between and 1.5 is randomly chosen 
for A; then, a number between and — v'Ada is chosen for 52- 



Br {h' , S' — >■ invis) = 

[T jh', S' AA)] [+ r (/^^ S' ^ AAAA)] 



(35) 



We calculate f ^ for each choice of the CxSM parameters 
to determine if the ATLAS experiment is sensitive to 
decays of the h' or S' eigenstates to dark matter. 



V. RESULTS 

In addition to fixing the Z2 breaking parameter ai = 
10~^, we also choose fixed representative values of some 
of the CxSM parameters for simplicity. We make the 
following choices for the parameters: 

• the dark matter mass is 10 or 100 GeV; 

• d2 {Mz) is fixed to 0.2, 0.5, or 0.9; 

• the singlet vev, x = 10, 100, or 1000 GeV. 

The values for the dark matter mass and the coupling 
were motivated by the study of the Z2a;SM in pS] 
which found a dark matter self-coupling of order 0.1- 
1.0 to be most interesting for satisfying vacuum stability 
while avoiding problems with non-perturbativity when 
the dark matter mass is between 10 and 100 GeV. The 
chosen order-of-magnitude values for the singlet vev re- 
sult in masses for the S' state which span a sufficiently 
large range that allows us to draw conclusions about the 
parameter space of the CxSM. We summarize in table [V| 
the values of the CxSM parameters chosen for our analy- 
sis. We will first present our results in detail for a single 
choice of parameters, and then present our general results 
for all those values in table IVl 



V.l. A Light Dark Matter Example 



We will use as our example the parameter set Ma = 
10 GeV, X = 100 GeV, and [Mz] = 0.2. The restric- 
tions on the CxSM parameters from the vacuum stability 
and RG analysis vary greatly with the choice of the cutoff 
scale A (described in section III I. This is demonstrated 



in the plots of (5| vs. \d2 in fig. [sj In the left column, 
§2 > 0, and on the right 62 < 0. The latter choice may 
accommodate a first order EWPT, as indicated by the 

work of rni. 



In each plot, the solid line indicates the tree level 
vacuum stability requirement (5| < Xd2. The plotted 
points correspond to values of the parameters that satisfy 
the RG running coupling constraints in eqn. (18) (gray 



points) or the traditional Landau gauge one-loop effec- 
tive potential vacuum stability requirement in eqn. ( |12[ ) 
(black points). Values of the cutoff scale A are taken 
to be 1 TeV (top row), 1000 TeV (middle row), or 
10^5 GeV ~ A/gut (bottom row). 

Fig. [3] evinces all of the generic features of the vacuum 
stability analysis discussed in section \IU\ Even with the 
most generous cutoff scale, A = 1 TeV, the allowed values 
of 62 for a given A and d2 do not extend up to the tree level 
bound because the RG evolution of the couplings breaks 
the condition of eqn. (18) that 62 (^i)^ < A (fj.) ^2 (m) at 
some /i < A. When 62 < (right column), this leads 
to a runaway direction in the Landau gauge effective po- 
tential: hence, the effective potential limits (black points) 
closely match the RG coupling limits (gray points). How- 
ever, if > (left column), the potential still appears 
stable (in the Landau gauge) even if the condition on the 
RG evolution of S2 is not satisfied and so the effective po- 
tential bound closely matches the tree level requirement. 

Furthermore, as the cutoff scale increases, small val- 
ues of A are forbidden because of the appearance of 
deep minima along the /i-axis of the potential (or, al- 
ternatively, A (^) < 0), as in the SM. Large values of A 
are also forbidden because RG evolution results in non- 
perturbative values for the quartic couplings in violation 
of eq. [19| 

In the discussion that follows, we consider the more 
conservative but gauge-independent bounds on the pa- 
rameter space from the RG evolution of the quartic cou- 
plings in eqn. ( 18 ), rather than the gauge-dependent lim- 



its from the effective potential stability requirement in 
eqn. (12). The plots for 62 > have similar features so 



we do not show them here; instead we focus on the 62 < 
scenario because of the impact on the electroweak phase 
transition, discussed in section [TV.ll 

The allowed masses of the h' , S' eigenstates are shown 
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FIG. 3. Plots of (5| (Mz) vs. A (Mz) d2 (Mz). For all plots, Ma = 10 GeV, x = 100 GeV, and d2 (Mz) = 0.2. The tree level 
vacuum stability requirement, S2 < Ad2, is indicated with the solid line. Gray points satisfy the constraints on the running 
couplings, eqn. (18 1, while black points satisfy the effective potential vacuum stability requirement, eqn. (121, in the Landau 



gauge. In the left column we take ^2 > while in the right column S2 < 0. The cutoff scale A is 1 TeV (top row), 1000 TeV 
(middle row), or 10^^ GeV (bottom row). 



in fig. |4] for the same set of parameters in fig. [Sj Ma = 
10 GeV, X = 100 GeV, da (Mz) = 0.2, and also A = 
1 TeV and S2 < 0. Fig. |4] shows the constraints on the 
masses from LEP searches, EWPO measurements, and 
the RG evolution of the quartic couplings in eqn. (18 1. 



Darker colored points result in a singlet relic density 
that is above the la WMAP bound, i.e., Q.Ah? > 0.118. 
Lighter colored points correspond to saturation or under- 
saturation of the relic density, O^ft,^ < 0.118. We note 



that the relic density is (under) saturated when there are 
s-channel resonances in the annihilation cross section (for 
2Ma — Ms') or the 4-point interaction dominates (for 
Ma > Ms'). 

One important feature of note is that increasing the 
Higgs-singlet coupling S2 decreases the mass of the lighter 
eigenstate — the singlet-like M5/ here — according to 
eqn. ^ and increases the mixing angle (see eqn. ([s])) 
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FIG. 4. Results of the scan for Ma = 10 GeV, x = 100 GeV, da (Mz) = 0.2, and A = 1 TeV with < shown in the Ms> 
vs. Mf^i plane. Dark colored points oversaturate the relic density, while light colored points (under)saturate. The top left plot 
imposes only the RG coupling limits, eqn. (18 1. RG coupling limits plus either the LEP constraints (top right) or the EWPO 
constraints (bottom left), and finally all three (bottom right), are also shown. 



when the other parameters (A, ^2, a;) are fixed. This effect 
is responsible for three features observed in fig. |4} 

1. The EWPO constraints (imposed in fig. [4] bottom 
left) favor a light scalar with SM-like couplings. As 
indicated by the slight slope on the right edge of the 
allowed region, larger h' masses are allowed as Ms' 
decreases due to the increased mixing that allows 
S' to offset the heavier h' . 

2. It is possible for the heavier Higgs-like eigenstate 
to avoid the 114 GeV bound from LEP (see fig. |4j 
top right) as Ms' decreases from a maximum of 
Ms' > 30 GeV due to an increased mixing angle 
and reduced h' coupling strength to SM particles. 

3. Also regarding the LEP constraints, a signifi- 
cant number of points corresponding to M5' < 
20 GeV — 2Ma are eliminated because the decay 
S' — >■ AA is no longer allowed, resulting in a light 
scalar 5" with SM-like branching fractions in viola- 
tion of the LEP constraints. Increasing the h' mass 
above 200 GeV decreases the mixing angle, so the 



lighter Algi masses are once again allowed by the 
LEP constraint despite the SM-like branching frac- 
tions of the S' . 

Fig. I5] and fig. [6] are the same as fig. |4] but with 
A ~ 10 GeV and 10^^ GeV, respectively. Increasing 
the cutoff scale forces 62 smaller — as was shown in fig. [3] 
— and results in relatively larger masses for both the 
scalar eigenstates. As these plots show, for a 1 TeV cut- 
off scale the effects of the LEP and EWPO limits are 
roughly equivalent. However, as A increases, the vacuum 
stability and perturbativity requirements reduce the al- 
lowed regions of parameters and masses; of the points 
that remain at these higher cutoff scales, a smaller num- 
ber satisfy the EWPO constraint than the LEP bounds. 
Thus the EWPO constraint is more significant than the 
LEP bounds at higher cutoff scales. Indeed, fig. [6] shows 
that there are scalar masses that satisfy the RG evolu- 



tion requirement in eqn. ( 18 ) up to the GUT scale, and 



the LEP bounds, but not the EWPO constraints. 

We now show the impact of the dark matter direct 
detection limits and the condition for invisibly decaying 
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FIG. 5. Same as fig.jijbut with A = 10*' GeV. 



Higgs searches at ATLAS. Fig. [7] displays those values of 
the Higgs-hke and singlet-Uke scalar masses that satisfy 
eqn. ( jlS] ) (top left — the same as the top left plot in 
fig.|4]), plus the XENONIOO direct detection cross section 
bound (top right) or the requirement for ATLAS invisible 
Higgs decay searches (bottom left), and all three together 
(bottom right). Fig. [s] and fig. [9] are similar to fig. [t] but 
with A = 10'' and 10^ GeV respectively. 

The choice of parameters here is such that nearly all 
the points satisfy the XENONIOO bound on the direct 
detection cross section. Consequently, very few points 
satisfy the CRESSTTI or DAMA regions (and a 10 GeV 
dark matter particle is incompatible with the result from 
CoGeNT presented in [39 ). Therefore, here and in 
what follows in later sections, we impose the XENONIOO 
bound as a more conservative upper bound on the direct 
detection cross section. 

More restrictive is the requirement for sensitivity to 
an invisibly decaying Higgs at ATLAS: this condition 
prefers a lighter Higgs-like eigenstate for which the total 
decay rate is smaller and hence the h' — > AA decay has 
a larger branching fraction. Though the light S' eigen- 
state has a large branching fraction to dark matter, the 
mixing angle is too small to give a greater than the 



requisite 60%. The RG evolution bounds and the ATLAS 
sensitivity become mutually exclusive at higher A. 



V.2. All Results for Light and Heavy Dark Matter 



We now consider all values of the couplings, singlet 
vev, and dark matter mass (light being 10 GeV, heavy 
being 100 GeV) listed in table |v| Fig. 



10 shows the ef- 



fect of varying the model parameters on the scalar mass 
eigenstates when Ma = 10 GeV; only the RG evolution 
constraints of eqn. ( jlS] ) have been imposed with a 1 TeV 
cutoff. The top center plot in fig. [T0| is identical to the 
top left plots in fig. |4] through fig. [9j 

Varying the singlet vev x between 10 and 1000 GeV 
clearly has a greater effect on the singlet-like eigenstate 
mass than varying the singlet quartic self-coupling d2 be- 
tween 0.2 and 0.9, as is expected from eqn. ([s]). Further- 
more, the S' eigenstate, when it is the lighter state (as 
in the left and middle columns) , has a maximum allowed 
mass: the smaller eigenvalue of the mass matrix has, for 
fixed d2,x, a maximum value of d^x^ jl (plus loop cor- 
rections) even as A increases. The Higgs-like eigenstate 
has a maximum value because we have limited our scan 
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FIG. 6. Same as fig. [i] but with A = lO^^^ GeV. 



of A. Finally, the general trend for the relic density is 
oversaturation when Mh' , Ms' > 2Ma — due to an off- 
resonance s-channel scalar exchange in the dark matter 
annihilation cross section — and undersaturation when 
2Ms' ~ Ma or Ms' < Ma- 



Fig. 11 and fig. 12 show the inclusion of LEP and 



EWPO constraints (right columns), direct detection and 
invisible search constraints (left columns), and higher 
cutoff scales of the effective theory for Ma = 10 GeV 
and 100 GeV respectively. Most of the discussion in sec- 
tion |V.1| generalizes to the other choices of the param- 
eters. We summarize the main results from these two 
figures as follows: 



• Avoiding oversaturation of the relic density requires 
at least one scalar eigenstate to be lighter than the 
dark matter except in the vicinity of a resonance in 
the annihilation cross section. Thus, if any of the 
dark matter direct detection experiments — in par- 
ticular, DAMA/LIBRA, CoGeNT, or CRESST-II 
— unambiguously detects lighter {Ma — 10 GeV) 
or heavier {Ma — 100 GeV) dark matter, it would 
be natural to consider the possibility of other light 
scalars weakly coupled to the SM. 



• The RG evolution and vacuum stability analysis 
requires heavier Higgs- and singlet-like eigenstates, 
i.e., less mixing, at larger cutoff scales (central 
columns). This is primarily to avoid the runaway 
direction in the potential corresponding to ^2 < 
and 62 {^if > A (^) ^2 {n). 



• the CxSM requires the existence of additional new 
physics below the GUT scale, regardless of whether 
the dark matter is lighter or heavier. If the S" eigen- 
state is lighter than the h' state, the RG evolution 
constraints and the LEP and EWPO limits (right 
columns) become mutually exclusive at high cutoff 
scales. If the S' eigenstate is heavier than the h' , it 
is possible to satisfy LEP and EWPO constraints 
with a lighter Higgs, but having a heavier S' results 
in oversaturation of the relic density. 



• The regions of parameter space most favorable for 
the invisible decay channel in the ATLAS detec- 
tor correspond to light (10 GeV) dark matter, a 
light singlet-like eigenstate, and a low (1 TeV) new 
physics scale. 
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FIG. 7. Results of the scan for Ma = 10 GeV, x = 100 GeV, da (Mz) = 0.2, and A = 1 TeV with < shown in the Ms> 
vs. Mft/ plane. Dark colored points oversaturate the relic density, while light colored points (under)saturate. The top left plot 
imposes only the RG coupling limits, eqn. ( 18 1. RG coupling limits plus either the XENONIOO bound (top right) or satisfaction 
of the requirement for ATLAS invisible Higgs decay searches (bottom left), and finally all three (bottom right), are also shown. 



V.3. Discussion 

Putting together the trends from all of the plots in 
fig. ID through fig. ^ we conclude that in order for the 
CxSM to be natural, the singlet vev cannot be very large 
(1000 GeV) and the scale of new physics has to be at most 
10^^ GeV. Moreover, if ATLAS observes the invisible de- 
cay mode of the Higgs, the new physics scale has to be 
in fact much smaller, at most 10 TeV, and the singlet 
self-coupling c?2 cannot be too large. To illustrate, for 
a; = 10 GeV, the maximum value of A for points that 
satisfy all constraints (LEP, EWPO, relic density, and 
is 10^^ GeV. However, as fig. [ill shows, 



direct detection) is 10 GeV. However, as fig. 
none of these points satisfy the requirements for the AT- 
LAS Higgs search via invisible decays. For x = 100 GeV, 
the maximum A for points that satisfy all constraints is 
about 100 TeV for = 0.2, 10 TeV for da = 0.5, and 
slightly greater than 1 TeV for da = 0.9. When we con- 
sider those points that also satisfy the requirements for 
the ATLAS Higgs search via invisible decays, the maxi- 
mum A is 10 TeV for da = 0.2, approximately 5 TeV for 
da = 0.5, and no points for da = 0.9. For x = 1000 GeV, 



there are no points that satisfy all constraints at any cut- 
off scale. 

It is interesting, then, to consider the additional im- 
pact of specifying the mass of the Higgs-like scalar to the 
range suggested by the recent ATLAS and CMS resultf]^ 
If the mass of the Higgs-like scalar is 125 GeV and AT- 
LAS does not detect the invisible decay mode, then the 
region near a; = 10 GeV is favored (the singlet-like state 
is very light) and the maximum A is '-^ 100 TeV|^ On 
the other hand, if the Higgs-like scalar mass is 125 GeV 
and ATLAS does detect the invisible decay mode, this 
seems to be compatible only with a singlet vev near 



^ We note that the significance of the reported excess will vary with 
the singlet-doublet mixing angle that affects the production cross 
section and the value of the dark matter mass that could allow 
additional decay channels to open. 

In these scenarios for which ATLAS is not sensitive to the in- 
visible decay mode, the Higgs branching fraction to dark matter 
may nonetheless be large and the branching fractions to visible 
final states reduced. This may be in conflict with the results 
from ATLAS and CMS, as pointed out in [9] for the ZaiSM. 
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FIG. 8. Same as fig. [t] but with A = 10*^ GeV. 



X = 100 GeV and d2 = 0.2. In this case, the singlet 
is again very light, with mass of order 20-25 GeV, and 
the CxSM requires new physics at 5 TeV, well within 
reach of the LHC. 



VI. CONCLUSIONS 



Two of the significant outstanding questions in par- 
ticle physics are the nature of dark matter and the en- 
ergy scale associated with physics beyond the Standard 
Model. Both of these topics are relevant for general sim- 
ple scalar extensions of the SM including the model we 
have studied, a complex scalar singlet extension. In the 
CxSM, through spontaneous and soft breaking of a global 
?7(1) symmetry, we obtain a massive stable dark matter 
candidate and two scalars which mix at the minimum 
of the potential, one Higgs-like and the other singlet- 
like. Rather than a traditional gauge-dependent vacuum 
stability analysis of the one-loop effective potential, we 
have chosen to place constraints on the parameters of the 
CxSM using a gauge-independent analysis of the renor- 
malization group evolution of the quartic couplings, mo- 
tivated by requiring the potential to be bounded below 



for the vacuum to be absolutely stable (metastable vacua 
may also be allowed, but we have not considered this 
possibility in our analysis) . Our analysis shows that con- 
straints on the RG running of the couplings gives results 
quite similar to the traditional vacuum stability analysis 
of the Landau gauge effective potential when the Higgs- 
dark matter coupling, 82, is negative as can be favorable 
to an EWPT; for ^2 > 0, the RG running constraints 
may be more conservative than is strictly necessary for 
vacuum stability. 

We have also considered constraints from relic density 
measurements, the electroweak phase transition, LEP, 
EWPO data, and dark matter direct detection experi- 
ments. Additionally, we have considered the sensitivity 
of the ATLAS experiment to a scalar which decays invisi- 
bly to dark matter. We have found that if the scale of new 
physics (the effective theory cutoff of the RG evolution 
of the CxSM parameters) is a TeV, then it is possible to 
satisfy all phenomenological constraints with a Higgs-like 
scalar mass in the region allowed by recent results from 
ATLAS and CMS and a light singlet-like scalar with a 
mass roughly twice the dark matter mass or less. Un- 
der these conditions, the mixing between the Higgs and 
the complex scalar singlet is very small, so the singlet- 
like scalar couples very weakly to SM particles. If the 



17 



Ma=10 GeV, x=100 GeV, ^2=0.2, A=10^15 GeV M^^IO GeV, x=100 GeV, d2=0.2, A=10^15 GeV 

30 ^ 



20 - 



15 




> 20- 



15 




80 100 120 140 160 180 200 220 80 100 120 140 160 180 200 220 

Mh' Mh' 
Ma=10 GeV, x=100 GeV, d2^0.2, A^IO'^IS GeV Ma=10 GeV, x=100 GeV, c/2=0.2, A=10^15 GeV 



30- 



20- 



15 



30- 



20- 



15 



80 



100 120 



140 160 



180 200 220 



80 100 120 



140 160 



180 200 220 



FIG. 9. Same as fig. [t] but with A = lO^^^ GeV. 



dark matter is light (10 GeV), the ATLAS detector may 
have sufficient sensitivity to a Higgs that decays to dark 
matter with a large branching fraction. 

If new physics does not appear at the TeV scale but 
instead arises at higher scales (10*^ GeV or the GUT scale, 
10^^ GeV), the CxSM is severely restricted by the com- 
bined vacuum stability and phenomenological bounds. In 
particular, limits from EWPO data are most in conflict 
with vacuum stability constraints (when the singlet-like 
eigenstate is lighter than the Higgs-like state) since the 
former generally favors a lighter Higgs while the latter 
requires a heavier Higgs at higher cutoff scales. 

With the continued operation of the LHC, a definitive 
statement on the existence of a SM-like Higgs and TeV 
scale new physics is expected. Determining whether or 
not the Higgs exists and whether or not it has SM pro- 
duction cross sections and branching fractions will shed 
light on the scalar sector of fundamental particle physics. 
Conclusive and consistent results from dark matter direct 
detection experiments will also provide information cru- 
cial for determining whether or not dark matter is a scalar 
particle and how dark matter couples to the SM. Until 
such results become available, our vacuum stability and 
phenomenology analysis has shown that further study of 



simple scalar extensions of the SM — 
CxSM with a light singlet-like scalar 



in particular the 
- is worthwhile. 
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Appendix A: One-Loop Potential & RGBs 



The one-loop potential is given by 



ViiK S,A) 



^E-^Tr{M;^(log 



Ml 



(Al) 

The sum i runs over the scalar, fermion, and vector bo- 
son contributions. The field-dependent mass matrices 
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FIG. 10. Here are shown the effects on the mass eigenvalues and the rehc density of changing d2 and x (with 82 < 0). Only the 
RG running constraints (A = 1 TeV) have been imposed. Left column (red): a; = 10 GeV; middle column (blue): x = 100 GeV; 
right column (green): x — 1000 GeV. Top row: = 0.2; middle row: d2 = 0.5; bottom row: d2 = 0.9. Dark colored points 
oversaturate the relic density, while light colored points (under)saturate. 



are given in eqns. (A2 1-(^11 ). The number of degrees of 
freedom associated with each contribution, n^, are given 
in eqn. (yll2|, and the numerical factors are given in 
eqn. (AlSl. In the scalar sector, we include the contribu- 
tions from the would-be Goldstone bosons. We use the 



notation g and g' for the SU{2)l and U{1)y gauge cou- 
plings, respectively, and yt is the top quark Yukawa cou- 
pling. Due to the smallness of the other fermion Yukawa 
couplings, we exclude them from the one-loop effective 
potential. 



Diag 



Mlih,S,A)., 



V 



mlhih,S,A) mls{h,S,A) ml^{h,S,A) 
mlg{h,S,A) mls{h,S,A) ml^{h,S,A) 
ml^{h,S,A) mlA{h,S,A) m\^{h,S,A) 



(A2) 
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FIG. 11. Here are shown the allowed mass eigenvalues (Ma = 10 GeV) after imposing the RG evolution requirement in 
eqn. ( |18[ ) (central column) plus the LEP and EWPO constraints (right column) or, alternately, the direct detection bound from 
XENONIOO and the invisible decay requirements (left column). The effect of increasing the cutoff scale is also seen: A = 1 TeV 
(top row), 10® GeV (middle row), and lO'^^ GeV (bottom row). Note that the top central plot corresponds to overlaying all 
plots in fig. |10[ The three colors indicate the value of x as in fig. |10[ The appearance of discrete bands for a given x is the 
result of the three different choices for d2. 



where 









3x3 


(A3) 


mlh {h, S, A) ^ 








(A4) 


mis {h,S,A) = 








(A5) 


^hA {h, S, A) = 
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- T 
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^AA {h, S, A) = 






ViA^) 


(A9) 
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2,2 



— -irh 

fcrmion ~ 2 



vector 



^scalar — 1 1 ^fc: 



= -2, riv 



= 3 



^scalar — 2 ' ^fermion — 2 ' '^vector — g 



(AlO) 

(All) 
(A12) 
(A13) 



The /3 and 7 functions are defined as 



fi d log Zy 
2 d^I 



(A14) 



for some coupling or mass parameter X and some field Y with wave function renormalization Zy- The /3 and 7 
functions for the CxSM are shown in eqns. (A15]-(A25). Since there are no interactions between scalars involving 
derivatives, at one-loop order there is no contribution to "fh, 7s : or 7^ from scalar loops. Thus Js — 7a — a-nd 7^ 
is unchanged from the SM result. 
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/3a = ^ Ux' + 51 ~ 36yt + \g'' + ^g^g'' + f / ) + 4A7. (A15) 



f^S2= (2rf2<52 + 2^2 ^ 352A) + 2^2 (7/. + 7s) (A16) 

^ {5dl + 25j) + Ad2js (A17) 



167r2 

/3,„2 = (62<52 + Sm^A) + 2m27;, (A18) 

/5''i= TTT^ (&ic«2) + 26i7s (A19) 

1 



Pb2 - 77^^ (262^2 + 2m2^2) + 2^275 (A20) 

Pa, = ai7s (A21) 

7s= 7A = (A23) 

(A24) 

Also included is the running of the vacuum energy, 

dn 1 fbj bl 



Appendix B: Oblique Parameters 



Here we present the oblique parameters S, T, and U, in terms of the gauge boson propagator corrections. We use 
the usual shorthand Cu, — cos 6'^ = M-[y /Mz, = sin( 



aT : 
aS - 
aU : 



1 



(47r)' 
9" /44 



(c^ [Hzz (0) + 2i^nz^ (0)] - n^w (0)) 



(47r)' V^H' 



Re \^zz (0) - (Af|) 
(0) - n^/vv (M^,)] 



[Hz^ (M|) - Hz^ (0)] 



o2 ^2 



n 



Z7 



(^) 



[^zz (Ml) 



n 

Ml 



(Ml) 



(0)] 



n 



(M|; 



(Bl) 



When calculating the differences AO between the CxSM and the SM values of 5, T, and [/, all terms except those 
dependent upon the scalar masses cancel. The relevant terms are given in eqn. (B2). 
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The functions Fq, Fi, and F2 (Afi2 = Fi — F2) are given by eqn. {B3 1 

»i 



{Mf,Ml,q^) = / dz In [(1 - z) Mf + zM^ - z {1 - z) 



(B3) 



M. Bento, O. Bertolami, and R. Rosenfeld, 
"Cosmological constraints on an invisibly decaying 
Higgs boson," Phys.Lett. B518 (20 01) 276- 281, 
arXiv:hep-ph/0103340 [hep-ph] 

M. Bento, O. Bertolami, R. Rosenfeld, and L. Teodoro 

"Selfinteracting dark matter and invisibly decaying 

Higgs," Phys.Rev. D62 (2000) 041302, 

arXiv : astro-ph/0003350 [astro-ph] ' 

J. McDonald, "Thermally generated gauge singlet 

scalars as selfinteracting dark matter," Phys.Rev. Lett 



[15] 



88 (2002) 091304 arXiv :hep-ph/0 106249 [hep-ph] , 



J. McDonald, "Gauge singlet scalars as cold dark 

matter," \Phys.Rev. D50 (1994) 3637-3649, 

arXivThep-ph/0702143 [HEP-PH] 

C. Burgess, M. Pospelov, and T. ter Veldhuis, "The 

Minimal model of nonbaryonic dark matter: A Singlet 

scalar," ^N ucl.Phys. B619 (20 01) 709-728, 

|arXiv : hep -ph/0011335 [hep-ph] [ 

X.-G. He, T. Li, X.-Q. Li, J. Tandean, and H.-C. Tsai, 

"Constraints on Scalar Dark Matter from Direct 

Experimental Searches," Phys .Rev. D79 (2009) 023521 

^rXi V : 08 1 1 . 0658 [hep-ph] 

Y. Mambrini, "Higgs searches and singlet scalar dark 
matter: when XENO NIOO meets the LHC," 
| arXiv: 1108.0671 [hep-ph] , * Temporary entry *. 
W.-L. Guo and Y.-L. Wu, "Exploring the singlet scalar 
dark matter from direct detections and neutrino signals 
via its annihilation in the Sun," [arXiv: 1103.5606 | 
[hep-ph] , * Temporary entry *. 

X.-G. He, B. Ren, and J. Tandean, "Hints of Standard 
Model Higgs Boson at the LHC and Light Dark Matter 
Searches," arXiv : 1112 . 6364 [hep-ph] 
M. Kadastik, K. Kannike, A. Racioppi, and M. Raidal, 
"Implications of the 125 GeV Higgs boson for scalar 
dark matter and for the CMSSM phenomenology," 
'arXiv: 1112.3647 [hep-ph] 

J. R. Espinosa, T. Konstandin, and F. Riva, "Strong 
Electroweak Phase Transitions in the Standard Model 
with a Singlet," arXiv: 1107.5441 [hep-ph] 
S. Profumo, M. J. Ramsey-Musolf, and G. Shaughnessy, 
"Singlet Higgs phenomenology and the electroweak 



phase transition," JHEP 0708 (2007) 010 
arXiv : 0705 . 2425 [hep-ph] 

V. Barger, D. J. Chung, A. J. Long, and L.-T. Wang, 
"Strongly First Order Phase Transitions Near an 
Enhanced Discrete Symmetry Point," arXi v: 1 112. 5460 1 
_[hep-ph] 

V. Barger, P. Langacker, M. McCaskey, M. J. 
Ramsey-Musolf, and G. Shaughnessy, "LHC 
Phenomenology of an Extended Standard Model with a 
Real Scalar Singlet," Ph ys.Rev. D77 (2008) 035005 , 
|arXiv : 0706 . 4311 [hep-ph] ( 



D. O'Connell, M. J. Ramsey-Musolf, and M. B. Wise, 
"Minimal Extension of the Standard Model Scalar 
Sector," Phys.Rev. D75 (2007) 037701 



|arXiv : hep-ph/06 lior4^ [hep-ph]| 

[16] V. Barger, M. McCaskey, and G. Shaughnessy, 
"Complex Scalar Dark Matter vis-a-vis CoGeNT, 
DA MA/LIBRA and XENONIOO," Phys.Rev. D8T 
I (2010) 035019, arXiv: 1005.3328 [hep-ph] 
[It] V. Barger, P. Langacker, M. McCaskey, 

M. Ramsey-Musolf, and G. Shaughnessy, "Complex 
Singlet Extension of the Standard Model ," \Phys.Rev. \ 



I D79 (2009) 015018, arXiv:0811.0393 [hep-ph] 
[TS] P. Fileviez Perez, H. H. Patel, M. Ramsey-Musolf, and 
K. Wang, "Triplet Scalars and Dark Matter at the 
LHC," Phys.Rev. D79 (2009) 055024, arXiv: 0811 . 3957] 
I [hep-phT] 

[19] J. R. Forshaw, D. Ross, and B. White, "Higgs mass 
bounds in a trip let model," JHEP 0110 (2001) 007, 
arXiv:hep-ph/0107232 [hep-ph] , 

[20] M.-C. Chen, S. Dawson, and T. Krupovnickas, "Higgs 
triplets and limits from precision measurements," 
Phys.Rev. D74 (2006) 035001 , arXiv:hep-ph/0604102~| 

I '[hep-phT[ " 

[21] M. Cirelli, A. Strumia, and M. Tamburini, "Cosmology 

and Astrophysics of Minimal Dark Matter," Nucl.Phys~\ 
' B787 (2007) 152-175, arXiv: 0706. 4071 [hep-ph]j 

[22] M. Sher, "Electroweak Higgs Potentials and Vacuum 
Stability," Phys. Rept. 179 (1989) 273-418 

[23] M. Lindner, "Implications of Triviality for the Standard 
Model," \Zeit. Phys. C3r (1986) 295 

[24] M. Lindner, M. Sher, and H. W. 



[25] 



M. Lindner, M. Sher, and H. W. Zaglauer, "Probing 
Vacuum Stability Bounds at the Fermilab Collider," 
Phys. Lett. B228 (1989) 139 
M. Sher, "Precise vacuum stability bound in the 



[26] 



standard model," \Phys. Lett. B317 (1993) 159-163 
arXiv : hep-ph/9307342" 

J. A. Casas, J. R. Espinosa, M. Quiros, and A. Riotto, 
"The Lightest Higgs boson mass in the minimal 
supersymmetric standard model," Nucl. Phys. B436 
T1995) 
A. 



3-29, arXiv:hep-ph/9407389,j 

[27] J. A. Casas, J. R. Espinosa, and M. Quiros, "Improved 
Higgs mass stability bound in the standard model and 
im plications for supersymmetry," Phys. Lett. B34 2 | 
(1995) 171-179, arXiv:hep-ph/9409458 
~A. 



] 



J. A. Casas, J. R. Espinosa, and M. Quiros, "Standard 
Model stability bounds for new physics within LHC 



reach," \Phys. Lett. B382 (199'6) 374-382 
larXiv : hep-ph/9603227" 

[29] M. Gonderinger, Y. Li, H. Patel, and M. J. 

Ramsey-Musolf, "Vacuum Stability, Pcrturbativity, and 
Scalar Singlet Dark Matter^" JHEP 01 (2010) 053 , 
|arXiv: 0910. 3167 [hep-phJI 



23 



[30] T. E. Clark, B. Liu, S. T. Love, and T. ter Veldhuis, 

"The Standard Model Higgs Boson-Infiaton and Dark 

Matter," arXiv : 0906 . 5595 [hep-ph] 
[31] R. N. Lerner and J. McDonald, "Gauge singlet scalar as 

inflaton and thermal relic dark matter," 

arXiv : 0909 . 0520 [hep-ph] 
[32] S. Caron and A. collaboration, "SUSY searches at 

ATLAS," arXivrilbe. lOOS Thep-ex]| * Temporary 

entry *. 

[33] A. Collaboration, "Search for neutral MSSM Higgs 
bosons decaying to r~ pairs in proton-proton 
colhsions at = 7 TeV with the ATLAS detector," 
arXiv: 1107.5003 [hep-ex] , * Temporary entry *. 

[34] CMS Collaboration Collaboration, S. Chatrchyan et 

al, "Inclusive search for squarks and gluino s in pp 

colhsions at ^ = 7 TeV" |arXiv : 1 107 ■ 1279 [hep-ex] 
* Temporary entry *. 

[35] CMS Collaboration Collaboration, S. Chatrchyan et 
al, "Search for supersymmetry in pp collisions at 
^/s = 7 TeV in events with a single lepton, jets, and 
missing transverse momentum," ,arXiv: 1107. 1870 | 
I [hep-ex] [ * Temporary entry *. 

[36] N. Nielsen, "On the Gauge Dependence of Spontaneous 
Symmetry Breaking in Gauge Theories," Nucl.Phys. ] 

yilOl (1975) 173, 
^.TLAS Collaboration Collaboration, "Combined 
Standard Model Higgs boson searches with up to 2.3 
fb-1 of pp collisions at sqrts=7 TeV at the LHC," Tech. 
Rep. ATLAS-CONF-201 1-157, CERN, Geneva, Nov, 
2011. 

[38] CMS Collaboration Collaboration, "Combined 
Standard Model Higgs boson searches with up to 2.3 
inverse femtobarns of pp collision data at sqrt(s)=7 TeV 
at the LHC,". 

[39] C. Aalseth, P. Barbeau, J. Colaresi, J. Collar, J. Leon, 

et al, "Search for an Annual Modulation in a P-type 

Point Contact Germanium Dark Matter Detector," 

arXiv: 1106.0650 [astro-ph.CO] , * Temporary entry 
* 

[40] D. Hooper and C. Kelso, "Impli cations of CoGeNT's 

New Results For Dark Matter," |arXiv: 1106 1066 ~| 

I [ hep-ph] , * Temporary entry *. 

[41] CoGeNT collaboration Collaboration, C. Aalseth et 
al, "Results from a Search for Light-Mass Dark Matter 
with a P-type Point Contact Germanium Detector," 
[Pftys.fe^.Lett. 106 (2011) 13130T||arXiv: 1002 .4703 | 
1^ [astro-ph.CO]" 

[42] R. Bernabei, P. BeUi, F. Cappella, R. CeruUi, C. Dai, et 

al, "New results from DAMA/LIBRA," Eur.Phys.J. \ 

I C 67 (2010) 39-49, arXiv: 1002. 1028 [astro-p h.GA] j "^ 

Temporary entry *. 
[43] C. Savage, G. Gelmini, P. Gondolo, and K. Freese, 

"Compatibility of DAMA/LIBRA dark matter 

detection with other searches," JCAP 0904 (2009) 010| 
arXiv: 0808. 3607 [astro-ph] * Brief entry^^'T 
[44] DAMA Collaboration Collaboration, R. Bernabei et 

al, "First results from DAMA/LIBR A and the 

combined results with BAMAfNal ," \Eur.Phys.J. C56 "| 

(2008) 333-355, arXiv: 0804. 2741 [astro-ph] ' 

[45] G. Angloher, M. Bauer, I. Bavykina, A. Bento, 
C. Bucci, et al, "Results from 730 kg days of the 

CRESST-II Dark Matter Search," arX iv: 1109.0702 | 

I [astro-ph.CO] * Temporary entry *. 



[46] B. Patt and F. Wilczek, "Higgs-field portal into hidden 
sectors," arXiv :hep-ph/0605188 [hep-ph] 

[47] Y. Zeldovich, I. Kobzarev, and L. Okun, "Cosmological 
Consequences of the Spontaneous Breakdown of 
Discrete Symmetry," Zh.Eksp.Teor.Fiz. 67 (1974) 3-11. 

[48] T. Kibble, "Topology of Cosmic Domains and Strings," 
J.Phys.A A9 (1976) 1387-1398 

[49] T. Kibble, "Some Implications of a Cosmological Phase 
Transition," Phys.Re pt. 67 (19 80) 183, 

[50] S. Profumo, L. Ubaldi, and C. Wainwright, "Singlet 
Scalar Dark Matter: monochromatic gamma rays and 
metastable vacua," Phys. Rev. D82 (2010) 123514, 
arXiv : 1009 . 5377 [hep-ph] , 

[51] H. H. Patel and M. J. Ramsey-Musolf, "Baryon 
Washout, Electroweak Phase Transition, and 



Perturbation Theory," 


JHEP 1107 (2011) 029", 


[arXiv: 1101. 4665 [hep-ph] 


* Temporary entry * 



[52] W. Loinaz and R. Willey, "Gauge dependence of lower 
bounds on the Higgs mass derived from electrow eak 
vacuum stability constraints," Phys. Rev. D56 (1997) 
I 741 6-7426, arXiv:hep-ph /9702321 [hep-ph] , 

[53] D. Boyanovsky, W. Loinaz, and R. Willey, "Gauge 

invariant Higgs mass bounds from the physical effective 
potential," Phys. Rev. D57 (1998) 100-111 
arXiv:hep-ph/9705340 [hep-ph] 

[54] G.-L. Lin and T.-K. Chyi, "Gauge independent effective 
potential and the Higgs mass bound ," \Phys.Rev. D6"0~| 
I (19 99) 016002, arXiv:hep-ph/9811213 [hep-ph]' 

[55] C. Wainwright, S. Profumo, and M. J. Ramsey-Musolf, 
"Gravity Waves from a Cosmological Phase Transition: 
Gauge Artifacts and Daisy Resummations," Phys. Re v. ~| 
I ' D8 4 (2011) 023521, arXiv: 1104.5487 [hep-ph], 

[56] T. Hambye and K. Riesselmann, "Matching conditions 
and Higgs mass upper bounds revisited," Phys. Rev. 
D55 (1997) 7255-7262, |arXiv:hep-ph/96 10272 
[he p-ph] I 

[57] K. Riesselmann and S. Willenbrock, "Ruling out a 

strongly interacting standard Higgs model," P hys. Rev. \ 

I D 55 (1997) 311-321| arXiv:hep-ph/9608280 [hep-ph] | 

[58] T. A. Chowdhury, M. Nemevsek, G. Senjanovic, and 
Y. Zhang, "Dark Matter as the Trigger of Strong 
Electroweak Phase Transition," ,arXiv: 1110.5334 | 

I [hep-ph] 

[59] Particle Data Group Collaboration, C. Amsler et al, 
"Review of Particle Physics," Phys. Lett. B667 (2008) ~| 
I 1-1 340 

[60] G. Belanger, F. Boudjema, A. Pukhov, and 

A. Semenov, "MicrOMEGAs 2.0: A Program to 
calculate the relic density of dark matter in a generic 
model," Comput.Phys. Commun. 1 76 (2007) 367-382 , 
|^Xiv:hep-ph/0607059 [hep-ph]] 

[61] XENONIOO Collaboration Collaboration, E. Aprile 
et al, "Dark Matter Results from 100 Live Days of 
XENONIOO Data," Phys. Rev. Lett. (2011) , 
arXiv: 1104.2549 [astro-ph.CO] , Long author list - 
awaiting processing. 

[62] XENONIOO Collaboration Collaboration, E. Aprile 
et al, "First Dark Matter Results from the XENONIOO 
Experiment," Phys. Rev. Lett. 1 05 (2010) 1313 02, 
arXiv: 1005.0380 [astro-ph. CD]] 

[63] CDMS-II Collaboration Collaboration, Z. Ahmed et 
al, "Results from a Low-Energy Analysis of the CDMS 

I^ Germanium Data," Phys. Re v. Lett.'i06j2011) \ 

I 131302, arXivilOll. 2482 [astro-ph. CO] | 



24 



[64] The CDMS-II Collaboration Collaboration, 

Z. Ahmed et al, "Dark Matter Search Resul ts from the 

CDMS II Experiment," Science 327 (2010) 1619-1621 

arXiv: 0912. 3592 [astro-ph.CO] 
[65] D. Hooper, J. Collar, J. Hall, and D. McKinsey, "A 

Consistent Dark Matter Interpretation For CoGeNT 

and DAMA/LIBRA," Ph ys.Rev. D82 (2010) 123509, 

arXiv : 1007 . 1005 [hep-ph]J 
[66] V. Barger, W.-Y. Keung, and G. Shaughnessy, "Spin 

Dependence of Dark Matter Scattering," Ph ys.Rev. 
I D 78 (2008) 056007, arXiv: 0806. 1962 [hep-ph] , 
[67] J. R. Ellis, A. Ferstl, and K. A. Olive, "Reevaluation of 

the elastic scattering of supersymmetric dark matter," 

Phys. Lett. B481 (2000) 304- 314, 

arXiv:hep-ph/0001005 [hep-ph] 
[68] G. Belanger, F. Boudjema, A. Pukhov, and A. Semenov, 

"Dark matter direct detection rate in a generic model 



wi th micrOMEGAs 2.2," Comput.Phys.Com mun. 180 | 
(20 09) 747-767, arXiv: 0803. 2360 [hep-ph] 

[69] A. Djouadi, J. Kalinowski, and M. Spha, "HDECAY: A 
Program for Higgs boson decays in the standard model 
and its supersymmetric extension," 
'T;om.put.Phys.Comm.un. 108 (1998) 56-74t 
arXiv : hep-ph/9704448 [hep-ph] 

[70] LEP Working Group for Higgs boson searches, 
ALEPH Collaboration, DELPHI Collaboration, 
L3 Collaboration, OPAL Collaboration 
Collaboration, R. Barate et al., "Search for the 
standard model Higgs boson at LEP," Phys.Lett. B565"| 
(2003) 61-75, arXiv:hep-ex/0306033 [hep-ex] 

[71] The ATLAS Collaboration Collaboration, G. Aad et 
al., "Expected Performance of the ATLAS Experiment - 
Detector, Trigger and Physics," ,arXiv: 0901 . 0512 | 
[hep-ex] 



